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ABSTRACT 



We report observations of a possible young transiting planet orbiting a previously 
known weak-lined T-Tauri star in the 7-10Myr-old Orion-OBla/25-Ori region. The 
candidate was found as part of the Palomar Transient Factory (PTF) Orion project. 
It has a photometric transit period of 0.448413 ± 0.000040 days, and appears in both 
2009 and 2010 PTF data. Follow-up low-precision radial velocity observations and 
adaptive-optics imaging suggest that the star is not an eclipsing binary, and that it is 
unlikely that a background source is blended with the target and mimicking the observed 
transit. Radial-velocity observations with the Hobby-Eberly and Keck telescopes yield a 
radial velocity that has the same period as the photometric event, but is offset in phase 
from the transit center by ~ —0.22 periods. The amplitude (half range) of the radial 
velocity variations is 2.4kms _1 and is comparable with the expected radial velocity 
amplitude that stellar spots could induce. The radial velocity curve is likely dominated 
by stellar spot modulation and provides an upper limit to the projected companion 
mass of Mpsini or b < 4.8 ± 1.2Mj up ; when combined with the orbital inclination, i or b, 
of the candidate planet from modeling of the transit lightcurve, we find an upper limit 
on the mass of the planetary candidate of M p < 5.5 ± 1.4Mj up . This limit implies that 
the planet is orbiting close to, if not inside, its Roche limiting orbital radius, so that it 
may be undergoing active mass loss and evaporating. 

Subject headings: Open clusters and associations: individual (25 Ori) - planets and 
satellites: detection - stars: individual (2MASS J05250755+0134243 / CVSO 30 / 
PTFO 8-8695 / PTF1 J052507.55+013424.3) - stars: pre-main sequence 



Introduction 



The Palomar Transient Factory (PTF) Orion project is a study within the broader PTF survey 
aimed at searching for photometric v ariability in the young Orion region, with the primary goal 
of finding young extrasolar planets (jvan Evken et al.l l201ll ). The project is based on a set of 
intensive high-cadence (~ 70-90 s) observations of a single 7.26 deg 2 field centered on the known 
25 Ori association, which lies in the Orion OBla region and has an estimated age of 7-10 Myr 
jBriceno et alJbood . l2007h . 



Typical young circumstellar disk lifetimes are on the order of 5-10 Myr (|Hillenbrandll2008l ) . and 
it is during this time that the bulk of the formation and migration of planets is e xpected to occur. 
The youngest exoplanets have been found via direct detection (e.g., LkCa 15, iKraus Ireland 
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2012 ; the free-floating planetary-mass object in the p Oph cloud, iMarsh et al.l |2010| ): however, 
little is known empirically about exoplanets during the first few millions of years of their lives, 
and the goal is to fill the observatio nal gap and be gin to provide constraints on theories of planet 
formation and evolution (see, e.g., I Armitagd 120091 ) . The stars in the 25 Ori/OBla region should 
be at or just beyond the point of di sk dissipation, cons istent with the high ratio of weak- lined to 
classical T-Tauri stars found there (IBriceho et al.l 120071 ) . We can therefore look for planet transits 
at the time when they may first become observable without their signatures being swamped by 
the extreme variability characteristic of the younger classical T-Tauri stars. In so doing, we can 
begin to investigate the frequency of planets at these ages; the timescales for their evolution; the 
timescales of their migration with respect to the star's evolution; and, through measurements of 
the transit depths, probe empirically their mean den sities and the extent of thei r atmospheres, 
which are expected to be inflated at these early stages (jFortnev &: Nettelmannll2010l . and references 
therein). 

Several surveys have been undertaken to search for close-in young exoplanets, though man 



are radial- velocity searches (e.g., Esposito et al. 20od: Paulson fc Yelda 2006 : Setiawan et al 



2008: Huerta et al.ll200S 



20071 : Miller et al.ll2008l : 



Crockett et al. 2011; 



Neuhauser et al 



2011 



Nguven et aDboid l. A transit seaich (e.g.. lAigrain et al 



many 
20071. 



has the advantage of being able to search many more 
stars simultaneously, down to fainter magnitudes and with a faster cadence. Combining transit pho- 
tometry with spectroscopic information can provide a wealth of information that is unavailable with 
non-transiting planets, from masses and radii to constraints on atmospheric composition. 



An outline of the PTF Orion project is given by Ivan Eyken et al.1 (120111 ). along with some of 
the first r esults conc e rning bina ry stars and T-Ta uri stars. The broader PTF survey is described in 
de tail by 



Law et al 



(|2009l ) and iRau et al.l (|2009l ). with a summary of the first year's performance 
by lLaw et al.l (|2010l ). Here we report a young planet candidate found orbiting a known M3 pre- 
main-sequence (PMS) weak-lined T-Tauri star (TTS) within the PTF Orion field. Our observations 
are outlined in Section [2 Section [3] discusses the photometric and spectroscopic results and some 
implications. The main conclusions are summarized in section [H 



2. Observations 



2.1. PTF Photometry 



PTF Orion data were obtained using the Palomar 48" Samuel Oschin telescope during the 
majority of the clear nights between 2009 December 1 and 2010 January 15, whenever the field was 
above an airmass of 2.0. All observations were in the R band, and of the 40 nights dedicated, 14 
yielded usable data, the remainder being lost primarily to poor weather. Light curves were obtained 
for ~ 110, 000 sources within the field, the top ~ 500 most variable of which were inspected visually. 
The observations and the PTF Orion d iffere ntial photometry pipeline and data reduction details 
are described fully in Ivan Eyken et al.l (|201ll ). Specifics of the more general PTF data reduction 
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Fig. 1. — Color-magnitude diagram following IBriceno et al.1 (Figure 7. 120051 ). highlighting PTFO 
8-8695 (CVSO 30) - the star symbol - in relation to the other T-Tauri stars discovered therein 
(see section [2TTj) . Open circles indicate WTTS; filled circles indicate CTTS. Solid lines indicate, 
from the top, 1, 3, 10 , and 30Myr isochrones, and the zero-age main sequence, according to the 
models of Si ess et al.l (120001 ) at a distance of 330 pc. Photometric measurements are reproduced 
from lBriceho et al.l (|2005l ). Reddening and extinction are neglected. 



pipeline are described by R. Laher et al. (2012, in preparation) and lOfek et al.l (|201ll ). Among 
the inspected sources, one, PTFO 8-8695 (2MASS J05250755+0134243), showed periodic shallow 
transit-like events with a shape and depth suggestive of a planetary companion, superposed on 
a larger-scale quasi-periodic variable light curve. The source has previousl y been identif ied as a 
we ak-line T-T auri star (WTTS) assoc iated with the Orion OBla region by lBriceho et al.l (CVSO 



30, 



20051 ) , and iHernandez et al.1 (|2007l ) . A summary of the main previously determined properties 



of the primary star is given in Table [H In Figure Q] we show a color-magnitude diagr am indicating 



PTF O 8-8695 in relation to the other T-Tauri stars in the OBla region discovered bv lBriceho et al 



(|2005l ). Although classified as a WTTS, it lies at the younger end of the distribution. This is 
consistent with its its relatively strong H a emission, which in fact place s it on the borderline with 
classical T-Tauri stars (CTTS) according to the classification scheme of IBriceno et al.l (|2005l ). 



The same field was again observed in the same way over seven clear nights between 2010 
December 8-17. The transit events were again evident with the same period and depth. Owing to 
improvements in the PTF image processing soft ware, better weather, and mitigation o f the CCD 



2011 



Ofek et al.ll20~ll\ the RMS 
noise floor for the brightest stars in the 2010 December data set improved from ~ 4 to ~ 3mmag, 



fogging effect seen in the previous year's data (see Ivan Evken et al. 
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Table 1. PTFO 8-8695 Stellar Properties 



Property 



Value 



Alternative designations 



a (J2000) 
5 (J2000) 
V 

2MASS J 
2MASS H 
2MASS K s 
Median R 
R range 
H a equiv. width 
Lil equiv. width 
Sp. Type 
T e s 
A v 

Luminosity 
Radius 
Mass (Baraffe/Siess) 
Age (Baraffe/Siess) 
Distance 



CVSO 30 

2MASS J05250755+0134243 

PTF1 J052507.55+013424.3 

05^25 m 07!55 

+01?34^24.3'/ 

16.26 mag a 

12.232 ±0.028 mag b 

11.559 ±0.026 mag b 

11.357 ± 0.021 mag b 

15.19 mag c 

0.17 mag (min to max) 

-11.40A a 

0.40 A a 

M3 (PMS weak-lined T-Tauri) a 

3470 K a 

0.12 mag a 

0.25 L a 

1.39i? a ' d 

O.44M /O.34M Q a ' e 

2.63Myr/2.68 Myr a ' e 

~ 330 pc (mean dist. to OBla/25 Ori assoc. ) a,f 



Briceho et all (|2005l ). 



ISkrutskie et al.1 (|2006|). 

c From PTF Orion data (this paper). 

d cf. smaller value implied by PTF Orion transit measurements - see Table [3j 



Reported bv lBriceno et al.l (|2005l ) based on comparison with lBaraffe et al.l (|1998l ) and lSiess et al 



(|2000h stellar models. 



Briceho et al.l (|2007l ) 
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with less evidence of systematic effects |j Figures 2 and 3 show the light curves obtained in the first 
and second years of observations with PTF. The quasi-periodic stellar variability is evident, along 
with sporadic low-level flaring, consistent with the star's young age and late spectral type. Visual 
inspection of the curves revealed the additional regular periodic transit-like signature, with a depth 
of ~ 3-4%, superposed on the stellar variability. We initially determined an approximate period 
(ks 0.45 days) and transit duration (~ 2 hours) by hand, in order to locate all the transit windows. A 
detailed discussion of the light curves and derived properties is given in Section [37TI The regularity 
and planet-like characteristics of the transit events led us to obtain follow-up observations. 



2.2. Keck Adaptive Optic Imaging 

The PTF imaging instrument is seeing limited and has a point-spread function (PSF) with 
typical full-width-half-maximum (FWHM) of 2.0", or 660 AU at the e stimated distance of the 



Orion OBla and 25 Ori associations (330 pc, iBriceho et al.l 120051 . 120071 ). We obtained adaptive 



optics (AO) images using the NIRC2 camera (PI K. Matthews) on the 10 m Keck II telescope in 
order to probe regions in the immediate vicinity of the star and to search for any sources (false 
positives) that could mimic the signal of a transiting planet, such as a nearby eclipsing binary. 
At R = 15.2 mag, PTFO 8-8695 is sufficiently bright for natural guide star observations and does 
not require the use of a laser for atmospheric compensation. We locked the AO system control 
loops onto the target using a frame rate of 41 Hz. The airmass was 1.29. We used the NIRC2 
narrow camera setting to provide fine spatial sampling (lOmaspix -1 ). Our observations consisted 
of 12 dithered H-band images (3 coadds per frame, 10 s per coadd), totaling 6 minutes of on-source 
integration time. Raw frames were processed by cleaning hot pixels, subtracting background noise 
from the sky and detector, and aligning and coadding the results. 

Figure H] shows the final reduced image using an effective field-of-view of 4.9" x 4.9", which 
corresponds to ~ 4.9 PTF pixels on a side. The image shows no contaminants, except for one faint 
source to the south east, at a separation of 1.8" (590 AU in the plane of the sky at the distance 
of Orion OBla/25 Ori). Our diffraction-limited images (FWHM ~ 80 mas) rule out additional 
off-axis sources down to a level of AH = 4.3, 6.4, 8.9, and 9.1 mag (3a) at angular separations of 
0.25", 0.5", 1.0", and 2.0" (83, 165, 330, and 660 pc) respectively. Assuming a color difference of 
Omag, a faint, blended, 100% eclipsing binary would have to be within AH ~ 3.5-3.8 mag of the 
primary star, or brighter, to mimic the observed transit depth at R of ~ 3-4%. At 6.96 mag (608 
times) fainter than our target, the one imaged contaminant is too faint to be a blended eclipsing 
binary capable of mimicking the observed transits unless it is extremely blue (R — H < —3.2); in 
that event, however, it would be unlikely to be stellar in origin in any case. Any other such binary 
would have to lie within 0.25" of our target in order not to have been detected. 



1 Independent differe ntial photometric analyse s of other PTF data, and respective precision levels, are discussed 
in lAgiieros et al.l 1)20111 ) and lLevitan et alj (1 20111 ). 
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Fig. 2. — Photometric light curve for 2009 December 1-2010 January 15 (Sections 12.11 and I3,l.ip , 
Gray regions indicate the transit windows, fixed at the measured transit period, width, and epoch 
of center-transit (To). Red points indicate data automatically flagged by the data reduction soft- 
ware as potentially non-optimal for various possible reasons (e.g., imperfect weather, evidence of 
contamination within the photometric aperture, etc.) 
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Fig. 3.— As for Figure for 2010 December 8-17. 
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Fig. 4. — Adaptive optics imaging of PTFO 8-8695, showing the full image field (left), and detail 
around the central star (right), at two different logarithmic brightness stretches. A very faint 
source is detected to the south east of the central star at a separation of 1.8", 6.96 mag fainter than 
the central source (indicated by the arrow, left panel); otherwise no contaminants are detected. 
(Section 

2.3. LCOGT, KPNO 4 m, and Palomar 200" Vetting 

Photometric observations were obtained in 2011 February with the Las Cumbres Observatory 
Global Telescope Network (LCOGT), using the 2m Faulkes North and South telescopes, and the 
0.8m Byrne Observatory telescope at Sedgwick Reserve (BOS). We were able to confirm complete 
transit detections on two separate nights with the BOS telescope with a clear filter (see Figure [5]). 
Observations in SDSS g' and i' with the three telescopes were inconclusive, and largely hampered 
by poor observing conditions and sparse phase coverage. The LCOGT observations enabled us to 
confirm the transit event and its period independently of the PTF data, and, in combination with 
the PTF data, we were able to establish a more accurate and up-to-date transit ephemeris. 

In addition we obtained low-precision followup radial velocity (RV) observations with the 
KPNO 4 m Mayall telescope and the Palomar 200" Hale telescope over 3 and 2 consecutive nights, 
respectively, to allow us to rule out a stellar-mass companion as the cause of the transits. The 
separate confirmation of the transit event and the lack of RV signature at the level of ~ lOkrns -1 
suggested we had indeed detected a sub-stellar object, and we pursued additional observations. 
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Fig. 5. — Transit observations obtained with the LCOGT BOS 0.8m telescope (clear filter), on 
the nights of 2011 February 8 and 9 (Section 12 .3|) . Gray regions again indicate the predicted 
transit windows using the period and To obtained with the PTF data. The second night shows a 
possible small flare at the end of the transit. Note that in addition to stellar variability, some of 
the systematic trends may also be air-mass related. 



2.4. HET and Keck Spectroscopy 

Following the AO vetting and low-precision RV followup, we obt ained Dopp ler radial ve- 
locity spectroscopy with both t he High Resolution Spectrograph (HRS; iTullI Il998l ) on the 9.2 m 
Hobby Eberl y Telescope (HET ; iRamsev et al.lll998l ). and the High Resolution Echelle Spectrom- 
eter (HIRES; IVogt et al.Ml994l ) on the 10 m Keck I telescope. The goal was to detect or place an 
upper limit on any signal due to reflex stellar motion caused by the orbiting companion. 



The queue-schedule operation mode of the HET (jShetrone et al.l 120071 ) allowed us to obtain 
four spectroscopic radial velocity observations very quickly after analysis of the low-precision RV 
and LCOGT photometric vetting data. The observations were timed on the basis of prior knowledge 
of the transit ephemeris to best constrain any orbital signature. 

Observations from the fiber-fed HRS spectrograph and calibrated using ThAr e xposures pro- 
vide radial velocity precision of better than 50 ms -1 over timescales of several weeks ( Bender et al 



2012 ). Seeing- induced PSF changes are minimal due to the image scrambling properties of optical 
fibers (Heacox 1987) and the HRS temperature is kept stable to ~ 0.01 C. We used the 15k reso- 
lution mode of HRS, with the red cross-disperser setting (316g7940), giving a wavelength coverage 
of 6114-9861 A. The use of this lower resolution setting enhances the efficiency of HRS without 
significantly degrading the information content in the stellar spectrum (since the target absorption 
lines were already known to be rotationally broadened). Each observation was 1200s in duration, 
except on UT 2011-02-21 which was 1600 s, and ThAr calibration frames were taken immediately 
after each to track instrument drift. Data for the four nights were reduced in IDL with a custom op- 
timal extraction pipeline that performs bias subtraction, flat fielding, order tracing and extraction, 
and wavelength calibration of the extracted one-dimensional spectra. A simultaneous sky-fiber was 
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available with this spectrograph setting, but we chose to avoid the additional complexity inherent 
in performing sky subtraction with a fiber-fed spectrograph, and instead masked areas affected by 
sky lines so that they were not used in calculating the RVs. The achieved signal-to-noise in each 
spectrum was approximately 40-50 per resolution element, peaking in the I-band. We identifie d 
the strongest absorption lines in an M2-type high resolution Kurucz stellar mode|§ ( Kurucz 2005 ). 
and used these to select lines for fitting in the measured spectra. 



For the Keck data, we observed PTFO 8-8695 using the standard procedures (e.g. jHoward et al 



2009) of the California Planet Search (CPS) for HIRES. The five observations over 10 days were 
each 500-600 seconds in duration and achieved a signal to noise ratio of 10-20 per pixel (SNR = 
20-40 per resolution element) in V-band. We used the "C2" decker (0.86" wide slit) for a spectral 
resolution of ~60,000. The total wavelength coverage was s=s 3650-7970 A. 

Owing to the faintness of the target and the relatively relaxed precision requirements, we 
employed in both cases thorium-argon (ThAr) emission lamps as the fiducial reference. This avoids 
the light-throughput penalty incurred by the higher precision common-path iodine gas absorption 
cell, which is better suited to brighter targets. 

Data reduction and analysis were performed independently at separate institutions for the two 
data sets, using independently written software, but following the same general principles. Visual 
inspection of the spectra readily revealed strong rotational line broadening ( Tasini* = 80.6 ± 
8.1kms _1 - see section [3.2.2|) . to the extent that only the very strongest of the photospheric lines 
were evident above the noise level. The low signal-to-noise (SN) and high stellar rotational velocity 
resulted in a very small number of lines from which to measure the radial velocity. Using a cross- 
correlation analysis on the HRS spectra to measure these broad, low SN lines yielded a velocity 
uncertainty of ~ 10 - 14kms _1 . Such an analysis is very sensitive to bad pixels or faint sky lines, 
which were difficult to remove in the low SN regime. In addition, a cross correlation strongly 
benefits by combining the information content of multiple lines with fixed relative spacings. In our 
data, viable lines tended to be spaced far apart: typically there were only one or two lines per 
spectral order, which negated much of the advantage of using a cross-correlation. Consequently, we 
adopted an alternative approach of fitting individual stellar lines with Gaussian functions, solving 
each line independently for its center, width, and depth. This fitting procedure was much less 
sensitive than cross-correlation to bad pixels or to the precise spectral window being analyzed. It 
also allowed us to mask out telluric absorption and strong, narrow night sky emission lines. 

The faintness of the target and the small semi-major axis of the putative companion's orbit 
meant that very high RV precision was neither anticipated nor required. Simple Gaussian+linear- 
term profiles provided an adequate fit (i.e., reduced \ 2 ~ 1) owing to the low SN and heavy 
broadening (and probable blending) of the lines, where more sophisticated models would have 
added unnecessary and unconstrained extra parameters. A weighted average of the shifts in the 



^http : //kurucz .harvard. edu 
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fit centroids with respect to the rest wavelengths of the stellar lines provided the required Doppler 
shifts. Since substantially differing SN in the various lines precluded a simple standard-deviation- 
based estimate of the measurement errors, errors were instead estimated by propagating the formal 
errors from the Gaussian fits to the lines. In the case of the Keck data, corrections were made for 
measured variation in the telluric lines to account for residual uncalibrated instrument drift, but 
these corrections were found to be negligible. A total of 14 absorption features spanning ~ 5300- 
7700 A were fit in the Keck data, and 14 spanning ~ 6100-8700 A in the HET data. The measured 
RVs are listed in table [5J 



3. Discussion 
3.1. Photometry 

3.1.1. Light curves and periodicities 

In order to model the transit events for a derivation of the transit and planetary candidate 
properties, the effects of the stellar variability in the light curves need to minimized - i.e., the light 
curves outside of transit needed to be whitened. After removing data points which are flagged or 
have large measurement errors, we fit a smooth cubic spline to all the PTF data which fall outside 
the transit windows (using the IDL imsLcssmooth functional), interpolating across the windows 
themselves. The fit is then subtracted (in magnitude space) from the entire light curve. Since 
most of the stellar variability occurs on longer timescales than the transits (with the exception 
of the occasional flares), this process yields the "whitened" light curve that has the majority of 
the stellar variability removed, leaving only the transits. Using a standalone version of the NASA 
Exoplanet Science Ins titute (NExScI) periodogram toolQ we calculate the Plavchan periodogram 



(jPlavchan et al.l 120081 ) of the combined whitened data sets, to provide a more accurate formal 



transit period measurement. We find a clear peak at 0.448413 ± 0.000040 d@ 

Figure shows the whitened data for the nights where any in-transit data were obtained, folded 
on the measured period. Short-term stellar variability is not corrected, and is probably the most 
likely explanation for those transits (particularly those from JDs 2455175, 2455192, and 2455545) 
which deviate significantly in shape from the general form of the others. This may be caused by 
a combination of low-level flaring, residual disk occultation, and the companion transiting small- 
scale stellar surface features. JD 2455192 appears to show a flare immediately after the transit, 



' 5 From the IDL Advanced Math & Stats module; see http://www.ittvis.com/idl 
4 See http : // exoplanetarchive . ipac . caltech . edu| 

5 We make the assumption that there is no phase shift in the transit timing between the first and second years, 
such that the datasets can be combined to give a year-long time baseline, and thus a very precise period estimate. If 
we disregard this assumption, the first year's data alone gives the best period estimate, P — 0.4486 ± 0.0010 d. 
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and also likely the tail of a second flare mid-transit, though the gaps in the data prevent a clear 
interpretation. JD 2455545 appears to show an early transit egress; however, the mid-transit 
variation and the disparity in comparison to other transits in the same year are suggestive of the 
above mentioned variability effects, which may confuse and mask the true egress time. 

The original "un- whitened" light curve is dominated by stellar variability (see Figures [2] and EJ. 
If we assume that the large-scale variability is caused primarily by spot modulation, and that the 
relative effect of the transit events is negligibl e, we can us e the un-white ned data to investigate the 



stellar rotation. Lomb-Scargle periodograms (|Lomblll976l ; IScargldll982l ) of the two years' data are 



shown in Figure[71 Where the Plavchan periodogram is designed for finding regular periodic features 
of arbitrary but constant shape (such as a transit), a Lomb-Scargle periodogram is better suited 
for finding periodic modes in a more complex signal such as the quasi-periodic stellar variation in 
our data: here, the signal appears not to repeat in an exact fashion and therefore does not fold 
well at any period (likely because of phase shifts due to changing spot features). 

A strong peak is found at 0.4481 ± 0.0022 d, in fact matching well with the transit period. 
Another peak is seen at 0.9985 ± 0.0061 d, corresponding closely to the sidereal or solar Earth day, 
and therefore most likely an artifact resulting from the observing cadence. The other peaks all 
appear to be aliases of these two periods. To confirm this, we modeled the data by creating an 
artificial light curve from two summed sinusoids with the same two periods, providing power at the 
frequencies of these two peaks. We superimposed artificial transits modeled as a simple inverted 
top-hat function, with the same depth, width, and ephemeris as the real transits. Allowing the 
phases and relative amplitude of the 0.4481 d and 0.9985 d signals to vary as free parameters, and 
requiring that the total amplitude of the light curve remain approximately the same as that of 
the actual data, we performed a least-squares fit to the real periodogram to assess how well its 
structure could be reproduced. Figure [7] shows the results, with a good match between model and 
data. Repeating the test with the 0.4481 d sinusoid omitted (i.e., summing only the one-day signal 
and the transits) gave a poor fit, implying that the form of the periodogram cannot be explained 
by the transits alone. Since there is clearly strong correlated out-of-transit variability associated 
with the star, and there are no other fundamental periods evident in the periodogram, the 0.4481 d 
signal appears to be the only likely period for the star. 

There are three other notable aliases of the presumed stellar rotation period, P*, at 0.3092 ± 
0.0010 d, 0.8126 ± 0.0070 d, and 4.43 ± 0.31 d (identified in the figure). The larger of these is 
substantially above the upper period limit implied by the measured f r sini* of the star (Section 
13X21) . P* < 27rP*/(v r sini*) = 0.671 ± 0.092 d. Though it is unlikely that the star would be 
coincidentally rotating at an alias of the transit period with the Earth's rotation, we repeated the 
modeling experiment with the model stellar rotation modified to match the r emaining two aliases 



to en sure that none could in fact be the true fundamental stellar period (see iDawson Fabrycky 



2010). Similar periodograms were obtained, with peaks at similar locations but with differing 
strength ratios; none were as good a fit to the periodograms based on the real data. We therefore 
conclude that the star is co-rotating or near co-rotating with the companion orbit. 
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Fig. 6. — Whitened light curves for nights where in-transit data were obtained, folded on the transit 
period (Section I3.1.11 left and right panels show first and second year's data respectively). Flux 
values are normalized to unity outside of eclipse, and nights are offset vertically in increments of 0.1 
for clarity. Transits are further distinguished in alternating black and gray, and crosses (red in the 
online color version) indicate data flagged as potentially compromised. The Julian day on which 
each transit occurred is indicated for comparison with Figures [2] and [3j The light gray regions 
indicate the transit windows. Note that partially covered transits which little or no data on one 
or other side of the transit window are likely to suffer from poor stellar-variability correction and 
show significant systematic error. 
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Fig. 7. — Lomb-Scargle periodogram of the non-whitened light curves from the two years of PTF 
Orion observations, compared to simple synthetic models (Section I3.1,ip . The presumed stellar 
rotation period is marked with a thick arrow (red in the online color version); a second distinct 
one-day period is marked with a thin (blue) arrow. The more prominent aliases of these periods are 
indicated with vertical lines: thick (red) marks indicate peaks which are predominantly produced 
by aliases of the presumed stellar rotation period; thin (blue) marks indicate peaks predominantly 
produced by aliases of the one-day period. Dotted (magenta) lines indicate peaks produced by a 
combination of the two. 
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We note that it is difficult to model the observed transit events with star spots alone, partic- 
ularly in the 2009 December/2010 January data: the short transit duration relative to the period, 
the flat bottom of the transits, and the sharp ingress and egress, are all more characteristic of a 
transiting object, whereas spots tend to cause smoother, more sinusoidal features as their projected 
area changes as they rotate across the stellar disk. It is more likely that we are seeing in the total 
lightcurve the effects of both star spots and a transiting object. 



3.1.2. Transit Fitting 

We fit transit models to the whitened light curves using the IDL-based Transit Analysis Package 



(TAP, iGazak et al.l l201ll ; Eastman et al, in preparation) . TAP employs Bayesian Markov Chain 



Monte Carlo (MCMC) tec hniques to explo r e the fitting parameter space based on the analytic 



transit light curve models of LM andcl Agoll (|2002l ). The package incorporates white- and red- noise 



parameterization ([Carter Sz Winnl |2009| ). allowing for robust estimates of parameter uncertainty 
distributions. We fit only the transits from nights where a complete transit was observed with 
adequate coverage both before and after the transit window, since on these nights the spline fit 
to remove the stellar variability is reasonably constrained on both sides of the transit; on nights 
where there is partial transit coverage, the spline fit tends to diverge during transit time. We also 
reject the transits mentioned in section I37L1 where stellar variability appears to strongly affect the 
shape of the transit, and in addition, the transit from JD 2455544, where the spline subtraction 
was particularly uncertain and rather sensitive to the tightness of the spline fit. Thus we fit four 
of the transits from the first year's data (JDs 2455201, 2455202, 2455205, and 2455211), and three 
from the second year (JDs 2455539, 2455540, and 2455543). The folded data are shown in Figure 

E 

For the transit fitting, we hold the orbital period, P, fixed throughout to the transit period 
determined above. The eccentricity, e, is fixed at zero, since it is difficult to constrain in the absence 
of a secondary eclipse. White and red noise levels are allowed to float separately for each individual 
transit, as are linear airmass trends. The remaining parameters are allowed to float, but are tied 
across all transits. They include: the epoch of transit center, To (we assume there are no transit 
timing drifts between the two years); the orbital inclination, i or b', a/R*, where a is the orbital semi- 
major axis, and R* is the radius of the primary; R p /R^, where R p is the companion radius; and the 
linear and quadratic limb-darkening coefficients for the primary star. In addition, we constrain the 
fitting such that R p /R* < 1, and white and red noise components are less than 4% (the depth of the 
transit). Parameters are calculated by fitting Gaussians to the dominant peak in the probability 
density distributions for each parameter resulting from the MCMC analysis. The center location 
of the Gaussian is taken as the parameter value, and measurement errors are estimated as the 
dispersion of the fit. We find the limb-darkening coefficients are largely unconstrained by the data. 
The fit is shown overlaid on the folded data in Figure [HJ the corresponding parameters are listed 
in Table El 
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Table 2. Differential Radial Velocity Measurements 



HJD 


RV 


J- V -J.Y.. kJ K..V.J I S \j 




kms -1 




2455613.668022 


1.81±0.64 


HET 


2455615.649694 


2.41±0.64 


HET 


2455616.640274 


-2.36±0.64 


HET 


2455623.622875 


0.55±0.64 


HET 


2455663.744361 


-1.53±0.91 


Keck 


2455670.747785 


-0.01±0.96 


Keck 


2455671.755651 


-1.38±1.03 


Keck 


2455672.757479 


0.10±1.02 


Keck 


2455673.770857 


2.83±1.33 


Keck 



Note. - RV uncertainties listed are formal uncertainties, and do not account for systematic 
effects. The offset between the two data sets is arbitrary, and set by shifting each set so that its 
mean is zero, excluding one outlier data point - see Sections 12.41 and 13.2.11 
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Fig. 8. — Folded light curve for the two years of observations combined, after removing stellar 
variability (Section 13. 1 . If) . The first year's data (2009 December - 2010 January) are plotted 
in black; the second year's (2010 December) are plotted in orange. A change in transit shape 
between the two years is evident. The best transit fit to both years combined is over-plotted on the 
assumption that R p < (Section l3.1.2p . Limb darkening is neglected, and error bars are omitted 
for clarity; the median photometric error is 0.0046, with an out-of-eclipse standard deviation of 
0.0056. 
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Table 3. System Parameters 



Parameter 



Value 



Measured 

P 

^orb 

a/R* 
Rp/R* 
T (HJD) 
v T sin 

Derived 



R* 
R p 

M p sin z orb 



0.448413 ± 0.000040 d 
61.8 ±3.7° 
1.685 ± 0.064 
0.1838 ± 0.0097 
2455543.9402 ± 0.0008 
80.6 ± 8.1 kms" 1 



0.00838 ± 0.00072 AU 
= 1.80 ± 0.15 _R a 
1.07 ± 0.10 Rq 
1.91 ± 0.21 R Jup 
< 4.8±1.2Mj„ b 



J up 



M p <5.5±1.4Mj 



up 



Note. — Summary of parameters determined in this paper for the PTFO 8-8695 system. Quan- 
tities are P - orbital period; i or b - orbital inclination; - inclination of stellar rotation axis; a - 
orbital semi-major axis; i?* - stellar radius; Rp - planet radius; To - epoch of transit center; v T - 
stellar equatorial rotational velocity; M p - planet mass. 



a From Kepler's third law, assuming stellar mass M, 
M p < M*. 

b 



0.39 ± 0.10 M Q (jBriceno et alJl2005l ). and 
Upper limit derived from measured RV semi-amplitude. 



-19- 



There is significant variation in the light curve from transit to transit, as can be seen in Figure 
[6j This can reasonably be attributed to the companion passing across varying surface features - 
cold or hot spots, or perhaps flares - on the stellar photosphere. Such variations are a result of the 
planet tracing the stellar surface brightness profile as it traverses the stellar disk, and cannot be 
removed by the whitening process, which is only sensitive to the integrated brightness of the disk. 
Since PTFO 8-8695 is expected to be active and spotted, such variation in the transits is further 
suggestive of a genuine transit, rather than a background blend. 

It can also clearly be seen from Figures [6] and [81 however, that there is an overall change in 
the transit shape between the two years' data sets. Explanations for this remain speculative. Re- 
running the fitting process and allowing a change in both i?* and R p between the two years yields a 
decrease in stellar radius of ~ 10% (2.1a) from one year to the next (with no measurable change in 
companion radius). Such a large change in the stellar radius in such a short time period is unlikely, 
and would likely manifest itself in observable rotation rate changes not seen in the data. A change in 
transit shape could also arise in principle from a change in orbital geometry. For example, the host 
star is rotating quickly enough to exhibit significant oblateness, which may induce precession of the 
orbital plane and, therefore, changes in i or b if the orbital and stellar rotation al axes are misaligned ^ 
An additional planet in a different orbit could also produce a similar effect ( Miralda-Escude 20021 ). 



It is difficult, however, to explain how an inclination change can yield a transit that becomes 
both more grazing (longer ingress/egress) and deeper at the same time. Another explanation may 
be variation in star spot coverage (or limb-darkening, though with heavy spot coverage, the two 
effects become somewhat confused). In addition to short- term spot variations, there may be surface 



features which survive for much longer periods (see, e.g.. IMahmud et al.ll201ll ). Given the activity 
of the star and the extreme proximity of the companion, such features may be compounded by 
magnetic or tidal star/planet interactions that could give rise to varying hot or cold spots near the 
sub-planetary point on the stellar photosphere. This could affect the apparent shape of the transit 
in a systematic way. The planet's apparent proximity to the tidal disruption limit (see section [373]) 
may also be a factor: a tidally distorted shape, or transient rings or a tidal tail of evaporating 
material could all yield unexpected and possibly varying transit shape (and also cause the slight 
transit asymmetry seen in the second year's data). For the sake of simplicity, and given the lack of 
data to disentangle all of these possibilities, we here adopt the single combined fit to both years' 
data sets, although with the caveat that the variability may cause some systematic error in the 
measurements. 



3 Given the small orbital period, this may occur on short timescales. Following Miralda-Escude ( 2002h . we estimate 



the order of magnitude of the oblateness-induced gravitational quadrupole moment of the host star, J2, by scaling 
to the pre-main-sequence from a solar value of ~ 10 -7 -10 -6 according to J2 oc 7?^/(M»P«). Assuming that the 
stellar obliquity is small, and noting that the orbital angular momentum and stellar rotational angular momentum 
are similar for M p ~ 5 Mj up (so that the inclinations of the stellar rotation and orbital planes with respect to the 
mean plane are similar), this leads to a precession period of the orbital nodes on the order of tens to hundreds of 
days. 
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We note that o ur fit yields a s omew hat smaller stellar radius, i?* « 1.07i? Q , than that pre- 
viously reported by iBriceho et al.l (120051) fl.39 Rp). A ssuming their estimate of T e g = 3470 K 
is correct, interpolating Siess models (jSiess et al.l |2000| ) with this updated radius gives a slightly 
older age estimate f or the star, ~ 3.7 Myr vs. 2.7 Myr. Given the distance uncertainties in the 
Briceho et al.l (|2005l ) results, however, and the possibility of uncertainty in T e g- arising from heavy 
spotting in the stellar photosphere, the two radius estimates are probably not inconsistent 



3.2. Spectroscopy 

3.2.1. Radial Velocities 

Since the RV analysis is differential in nature, the RV offset between the HET and Keck 
data sets is arbitrary. We place them on the same approximate scale by shifting the RV zero 
points to the mean of the data for each data set. There are too few data points to create a 
periodogram to measure any periodicities in the data, but we can look for consistency with the 
previously determined transit period by folding the data on that period and looking for a coherent 
alignment of the data points. The result is shown in Figure O indeed, the data appear to phase 
well, showing a smooth and apparently sinusoidal variation, with the exception of one outlying 
data point at orbital phase <j> = —0.082 (Heliocentric Julian date (HJD) 2455615.64969), which is 
discussed further below. This outlier is excluded in calculating the offset between the data sets. 

To constrain t he mass of the companio n, we fit a Keplerian orbit model to the data using the 
RVLIN package bv lWright fc Howard! (|2009l ). The model includes 6 parameters: the period, P; the 
RV semi-amplitude, K; the eccentricity, e; the argument of periastron, uj; the time of periastron 
passage, T p ; and the systemic velocity, 7. We fixed P and T p to the values measured from the 
transit photometry. Since the \ 2 surface has multiple minima, we explore the fitting parameter 
space by running 10000 trial Keplerian fits (neglecting the apparent outlier point), with initial 
parameter estimates drawn at random from reasonable starting distributions, and selecting the fit 
with the lowest reduced x 2 (x 2 e d)- 

We find that the constraints from the transit ephemeris make it difficult to obtain a reasonable 
Keplerian fit. A good circular-orbit fit, with the eccentricity, e, fixed to zero, is prohibited by the 
constraint on To provided by the photometry. Such a fit must cross its central velocity at phases 
and 0.5, with minima and maxima at 0.25 and 0.75 (the quadrature points), and as a result appears 
significantly out of phase with the data (x 2 e d = 4-0; dashed line, Figured]). 



An eccentric fit goes some way to resolving the mismatch (x 2 e d = 1-1)> but yields an eccentricity 



7 Our ra dius estimate also de pends on our estimate of a, which depends in turn on the assumption that the mass 
estimate of Briceno et al. ( 20051 ) is correct, so the argument is somewhat circular, a, however, is relatively insensitive 



to stellar mass (a oc M 1 ^ 3 ), and the Siess models predict a negligible difference in mass at our older age. (In fact, at 
the estimated T e g, they predict a mass range smaller than the errors over an age range of 1-10 Myr.) 
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Fig. 9. — Differential radial velocity measurements obtained with Keck/HIRES and HET/HRS 
(Section I3.2.ip . Zero phase is chosen to correspond to the photometric center-of-transit time. The 
offset between the two datasets is chosen to match the mean RV of each. The lines indicate best 
Keplerian fits (excluding the outlier): dashed line - circular orbit, transit-center time, To, fixed 
to photometry (x^ed = 4.0); dotted line - eccentric orbit, transit-center time fixed to photometry 
(Xred = 1'1)> son d nne ~ a sinusoidal fit (equivalent to a circular orbit) at the same period, with 
phase free to float (Xred = 0-42). The eccentric fit is better than the fixed- To fit, but brings the 
companion to the surface of the star at periastron. The floating-phase sinusoidal fit gives the best 
result, suggesting that star spots either modify or dominate the Doppler RV signal. The outlier 
point may represent Rossiter McLaughlin effect due to the transiting companion. 
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of e ~ 0.492 (dotted line, Figure [9]). This places the fractional periastron distance at r per i/a = 
1 — e ~ 0.51, bringing the companion to the point of contact with the surface of the star as 
determined from the photometric transit fits (R*/a ~ 0.51). Such a solution may be improbable, 
given that in section [3T31 we note that the orbital semi-major axis appears to be close to the Roche 
limit. 



A better fit is in fact obtained by fitting a circular orbit model and removing the constraint 
on To, allowing the phase to float. This fit is over-plotted in Figure [U showing the fit is good 
0.42), but offset in phase from the transit ephemeris by —0.22 ± 0.04 periods. 



(^red 



We favor this better fit to the RV signal, which most likely arises because of spot effects 
modulated by the stellar rotation, where the amplitude of the spot effect is at least comparable 
to - if not much greater than - any true reflex Doppler signal from the companion. Similar spot- 
induced RV amplitudes for have already been observed fo r other T-Tauri stars in the optical (e.g., 
Mahmud et al.ll201ll ; lHuerta et al.ll2008l : iPrato et al.ll2008l ). Since the star appears to be co-rotating 
with the planet orbit, the period of such a spot-dominated RV signal would match the orbital period, 
but the phase would be arbitrary, depending on the longitudinal spot placement. Alternatively, 
spot and companion RV signals may both be significant and the two effects may compete (see 
Section [3J3]). 



Regardless of the astrophysical cause behind the measured RV signal, we can place an upper 
limit on the companion mass, and because of the good orbital phase coverage we can be reasonably 
assured that possible aliasing effects are not a concern. We assume that any companion-induced 
Doppler motion must be smal ler than the amplitu de of the measured variations, that M p <C M*, 
and that M* = 0.39 ±0.10M^ (IBriceho et alJbood ). and hence estimate M p sinz orb < 4.8±1.2M Jup 
(see, e.g., formula 3. lGaudi Winnll2007l ) . Taking our derived value of i or b we can directly estimate 
M p < 5.5 ± 1.4Mj up , comfortably within the planetary mass regime. The semi-amplitude (half 
peak-to-peak) of the combined measured RV variations is ~ 2.4 km s . This is separately confirmed 
by the two individual datasets, which both have good phase coverage and show similar amplitudes 
independent of the RV offset (as can be seen in Figure [9]). By comparison, since Doppler- induced 
RV semi-am plitude scales dir e ctly w ith M p , a 25 Mj up object on the planet/brown-dwarf boundary 
suggested bv lSchneider et al.l (|201ll ) would induce an llkms -1 signal; an object on the deuterium- 
burning limit (« 13Mj up ) would induce a 5.7km s _1 signal. It is unlikely, therefore, that the object 
is of more than planetary mass, and still less in the stellar mass ranged Since the spot distribution 
is likely to change with time, further observations to look for changes in the phase of the RV signal 
with respect to the transit ephemeris may provide more insight into the exact nature of the RV 
signal. RV observations in the infra red are also likely to suffer less from spot-induced noise, and 
so could also provide valuable further information. 



8 We do note, however, that there is a possibility that if the reflex Doppler signature and the spot signature 
are of comparable amplitudes, matching periods, and opposed in phase, they may destructively interfere, giving an 
artificially low RV amplitude. 
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We were unable to find any abnormalities in the observations regarding the outlying RV data 
point, and, thus, have no reason to reject it as a bad measurement. However, falling at an orbital 
phase <f> = —0.082, it lies within the transit window (c.f. Figure [8|), and its ano malous velocity 
could reasonably be explained as arising from the Rossiter-McLaughlin (RM) effect (|Rossiterlll924l ; 



McLaughlinlll924l ) . The RM effect would appear regardless of whether the main RV signal is caused 



by spot modulation or stellar reflex motion, since it is a result of an asymmetric distortion of the 
stellar absorption lines as an obscuration transits the unresolved stellar disk, blocking off regions 
with successively different rotational red-shifts. Indeed, given the rapid stellar ro tation, the RM 
effect should be expected to appear during the transit window. From equation 6 of lGaudi &: Winn 
(|2007l ) we can estimate the expected maximum amplitude of the effect, Kr, given v r sinz*, R p , and 
i?*. We find Kr ~ 3kms _1 , in good agreement with the offset of the outlier. The RM maximum 
also lies typically around ~ 1 /5— 1 /3 of the way between transit-ingress and transit-center for gas- 
giant planets, which is coincidentally where the outlier point is located. The sign of the offset of the 
outlier data point, falling prior to the transit center time, is in agreement with a prograde orbit, as 
we would expect if the star is in synchronous (or quasi-synchronous) rotation. It should be noted 
that the RV data point at <f> = +0.039 (HJD 2455671.755651) also lies within the transit window 
but, in contrast, does not appear to be an outlier. Lying closer to the transit center time, one 
would expect the RV offset to be smaller here; however, the exact form of the RM effect depends 
on the precise orbital geometry, and may be partially masked by the choice of offset between the 
two data sets. Further complication may also arise from the companion transiting photospheric 
surface features. Dedicated RV measurements during transit could provide valuable confirmation of 
the RM effect hinted at by the data, and help independently confirm the validity of the transiting 
planet candidate. 



3.2.2. Stellar Rotation 

The stellar rotation velocity, v r , measured from the spectroscopy provides an independent 
consistency check on the stellar rotation period. In order to estimate the projected rotation velocity, 
■u r sini*, we degrade the Kurucz synthetic model spectrum (see Section I2~4"l) to the resoluti on and 



sampling of the HRS, and rotationally broaden it using a non-linear limb-darkening model (| Claret 



2000 



Gravlll992l ) over the range of v r sin i* = 10 - 150 km s _1 . We then Doppler shift the broadened 
models, line by line, by cross-correlating against lines in our target spectrum. Subtracting a Doppler 
shifted model from the observed spectrum gives an RMS residual for that rotational velocity. 
Minimizing this residual gives the optimal rotational velocity measured for an individual line. We 
applied this analysis to 12 lines that were sufficiently deep to be visible above the noise for our full 
range of u r sini* models. Weak lines could not be used because at large u r sini* the line profiles 
became buried in the noise. The weighted average of v r sin measured for these lines, across all 
four HRS observations, was 80.6 ± 8.1 km s -1 . 



The assumed rotational period of the star and the radius measured from photometric transit 
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fitting (neglecting oblateness effects) independently imply an expected value for the unprotected 
rotational velocity of v r = 27ri?*/P* « 120 ± llkms -1 . Taken together, the measured v r sini* and 
the photometrically derived v r give an estimated inclination of the stellar rotation axis, i* = 42 ±7°. 
Comparing this with the measured orbital inclination, i or b = 61.8 ± 3.7°, we see that there is weak 
evidence for a possible misalignment between the orbital and stellar rotation axes (consistent with 
the possibility of detecting changes in the orbital inclination on relatively short timescales, as 
mentioned in Section f3.1.2j) . 



3.3. Implications 



In Figure [10] we show the radius and mass of the companion in relation to the currently known 
transiting exoplanets listed in the NASA Exoplanet Archive^ where we have indicated the mass 
of the candidate as an upper limit only. It clearly lies at the upper end of the gas-giant radius 
distribution, although not unprecedentedly so. A significantly inflated atmosphere is to be expected 
owing both to the system's very young age, and to stellar irradiation at the companion's exception- 
ally small orbital radius. Planetary atmosphere models are currently not wel l constrained i n thi s 
regi me due to the l ack of known exoplanets that are both young and close-in. Mar ley et al. (j2007i ) 
and iFortney et al.1 (120071 ) caution that evolutionary models should be treated with care at ages 
up to 10 Myr or more, being highly d e pend ent on initial conditions and the formation mechanism 
assumed. Indeed, ISpiegel Burrows! (|2012l ) suggest that comparison of atmospheric models with 
observations of exoplanets at s uch young ages may be a good way of distinguishin g between differ- 
ent formation scenarios. Both ISpiegel Burrows! (J2012J) and iMarlev et al.l (|2007l ) predict radii of 
around 1.6i?j up for a 5Mj up planet at 3 Myr, for non-irradiated 'hot-start' (gravitational collapse) 
model atmospheres; our error bars place our measured radius (1.9 ± 0.2 Rq) just above this value. 
Our radius lies substantially further above the post-formation cold-start (core collapse) isochrones, 
though comparison with the latter models is confused by uncertainty in the formation timescale, 
which is likely comparable to the age of the star. Accounting for stellar irradiation could increase 
the theoretical radius further. Inc reases of ^ 10% or more are typical at ~ 0.02 AU from a so- 



lar like star (Baraffe et al 



2010 



Chabrier et al.l 120041 ) . where the incident flux is similar to that 



on our planet candidate (which is closer in, but orbiting a lower-mass star). The companion's 
proximity to the Roche limit may also weaken the gravitational binding of the object, further in- 
creasing its expected radius. Various other explanations are also proposed to explain the excess 



'radius anomaly" th at is increasingly found in other gas-giant exoplanets (e.g.. iBaraffe et al. 



Chabrier et al. 



2010 



2011 



). Our estimated radius is therefore not unreasonable, but due to both the 



model and observation uncertainties, a robust comparison with theory is difficult. 

The apparent companion to PTFO 8-8695 also orbits close enough to its parent star that the 
consequent small size of its Roche lobe may be relevant. Below a certain mass threshold, it may 



: //exoplanetarchive . ipac . caltech.edu 
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0.1 1.0 10.0 

Mass (M Jup ) 

Fig. 10. — Radius vs. mass for the known transiting planets (Section 13. 3t data taken from the 
NASA Exoplanet Archive). PTFO 8-8695b is marked as an upper mass limit, highlighted by the 
large square. Iso-density contours are marked at 0.5, 1, 1.33 (Jupiter density), 5, and 10gcm~ 3 . 
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not have sufficien t self-g ravity to hold i tself t ogether, and thus may begin to lose mass. Following 
Faber et al.l (|2005l ) and iFord k, Rasid (|2006l ). the Roche radius, -RR oc h e > according to iPaczviiski 
(|197ll ). is given by 

^Roche = 0.462(M p /M,) 1 / 3 a < 0.462(M p , max /M,) 1 / 3 a « 1.92 ±0.16, (1) 

where M Pjmax is our previously estimated maximum mass for the companion. Comparing with the 
measured companion radius gives -R p /i?R chc It 0.994 ± 0.094. Framing the argument another way, 
we can rewrite the Roche formula to estimate the Roche limiting orbital radius, i?R oc he! i n terms 
of the measured companion radius to find a/aR oc h e % 1.008 ± 0.095 at the most - consistent with 
being at or within the Roche limit, within the errors. The planet may be sufficiently inflated that 
it fills its Roche lobe, and consequently may have lost mass in the past, or be in the process of 
losing mass (thus maintaining itself at the Roche limit). 



Conclusions 



We have detected transits from a candidate young planet orbiting a previously-identified co- 
rotating or near-co-rotating ~ 2.7Myr-old M3 weak-line T-Tauri star. Although we cannot com- 
pletely rule out a false positive due to source blends, we are able to rule out confusion at the level 
of AH w 4.3 mag beyond a separation of 0.25", and argue qualitatively that a false detection due 
to a blended eclipsing binary is unlikely. The companion is in an exceptionally ra pid 0.448413 d 
orbit , placing it among the shortest of the currently known exoplanet p eriods (cf. iDemory et al. 
201ll ; ICharpinet et alj|201ll ; iMuirhead et al.l 120121 ; iRappaport et al.ll2012l ). With an orbital radius 
only around twice the stellar radius, it appears to be at or within its Roche limiting orbit, with 
a/aR oc he < 1.01±0.10, raising the possibility of past or ongoing evaporation and mass loss. Perhaps 
the companion has been migrating, and losing any mass beyond its Roche lobe as it does so; or 
perhaps it is continually being inflated to fill its Roche lobe, with any material which overflows 
being stripped away. 

Although the transit photometry and the RV data both phase-fold on the same periods, there 
is an apparent offset in phase between them. The most likely explanation is that the RV signal is 
shifted or dominated by the effect of star spots; we therefore suggest an upper limit on the (incli- 
nation independent) companion mass of ~ 5Mj up based on the amplitude of the RV modulation. 
If it can be assumed that the object has had time to reach a stable (or quasi-stable) state - i.e., 
that mass loss rates are not too rapid, and that there have been no recent dramatic changes in the 
orbital geometry - then Roche limit considerations would imply a lower limit of a similar order, 
since anything much less would be unable to gravitationally bind the material within the measured 
companion radius. 

The data are complex enough that we cannot yet be certain of a planet detection. In favor of 
the planetary interpretation, however, we note that: 1) the photometric transit shape appears to be 
flat bottomed with sharp ingress and egress slopes, and is difficult to explain with a star-spot model 
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alone; 2) the transits appear highly consistent and periodic over a period of > 1 yr, which would be 
unusual for spots; 3) the RV signal places an upper mass limit well within the planetary regime; 4) 
a false positive due to a faint blended eclipsing binary is argued against by two observations which 
are suggestive of a transiting object associated with the primary T-Tauri star observed: the stellar 
rotation rate appears very close to the transit period, suggesting co-rotation; and the photometric 
variation from transit to transit is consistent with the notion of a planet transiting a heavily spotted 
T-Tauri star. 

Further spectroscopic and photometric observations are needed to provide a clearer picture. 
Since the star is such a fast rotator, the expected RM effect (hinted at in the data) may provide a 
valuable opportunity for confirmation of the candidate. Most RV noise sources will be constant in 
the timescale of one transit, so further RV observations during the transit window could potentially 
provide full sampling of the effect. This would help confirm the planetary interpretation of the data, 
and could provide further useful information on the system geometry. Revisiting the target to make 
further measurements of the RV phase with respect to the transit ephemeris would help confirm the 
spot-interference interpretation: if spots are significant, the phase offset is likely to change with time 
as the spots change. Infra-red RV observations, which are less sensitive to spot-induced noise, would 
also provide valuable further information, perhaps allowing for an unambiguous determination of 
the companion mass. 

Given the young age of the system, the planet candidate is likely hot, and so a secondary 
transit may also be detectable with more precise photometry, allowing constraints to be placed on 
the companion temperature. Finally, simultaneous multi-band photometry and spectroscopy (par- 
ticularly in the infra-red, where spot-effects should be lessened) could further help disentangle the 
signatures of star spots and companion. This would also help to establish the cause of the apparent 
overall change in transit shape between the two years' observing runs, which could result either 
from changes in the spot distribution or possibly, given the exceptionally short orbital timescale, 
changes in the orbital geometry. 

If our interpretation of the data is correct, the putative planet's youth and its unique proximity 
to its host star will make it a valuable object for helping inform our understanding of exoplanet 
formation. Its young age could have important i mplications for the mechanism by which it formed: 
conventional core-accretion formation models (jPollack et al.l 1 19961 ) occur on timescales of ~ 1— 
lOMyr, comparable wit h or longer than the age of the system; the much more rapid gravitational- 
instability mechanism (jBossI 119971. 120001) occurs on timescales orders of magnitudes shorter, and 
may thus perhaps be favored (see lBaraffe et alj|2010i . for a brief overview and comparison of the two 
mechanisms) . The comp anion's inflated atmospher e appears indicative of the 'hot start' models of 
Mar ley et al.l (|2007l ) and ISpiegel Burrows! (120121 ) , which are associated with gravitational insta- 
bility. At the same time, we cannot rule out that formation may yet be incomplete given the young 
age of the system. Neither can we rule out that the companion may be on the verge of evaporation 
and may not survive the lifetime of the star. Indeed, very few large planets are known to orbit 
small stars, and this object may well be transient. Much ambiguity remains at these young ages, 
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and more observation and analysis of the system is needed. If the planetary nature of the proposed 
companion is confirmed, the system could provide a wealth of new information for constraining 
dynamical and atmospheric evolution models for exoplanets. We therefore propose PTFO 8-8695 
as an object worthy of further careful study by the community. 
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